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We develop the quark coalescence model by including the quark number conservation in order 
to describe the hadronization of the bulk of the quark-gluon plasma. The scalings between hadron 
and quark phase space distributions are shown to depend on the transverse momentum. For hard 
quarks, our general scalings reproduce the usual quadratic scaling for mesons and the cubic scaling 
for baryons. For softer quarks, the inclusion of the quark number conservation leads to a linear 
scaling for the hadron species that dominates the quark number of a flavor, while the scalings of 
non-dominant hadrons are sensitive to the coalescence dynamics. For charm mesons, we find that the 
distribution of soft D mesons does not depend on the light quark distribution but the distribution 
of soft J/ip mesons is inversely correlated to light quarks. 

PACS numbers: 12.38.Mh, 25.75.Ld, 25.75.Nq 
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1 Introduction 

For the dense matter created in ultra-relativistic heavy 
ion collisions such as those at the Relativistic Heavy Ion 
Collider (RHIC) and the Large Hadron Collider (LHC), 
the quark coalescence model provides a hadronization 
model. When the quark-gluon plasma reaches the 
hadronization hypersurface after expansion, a quark and 
an anti-quark can form a meson while three quarks can 
form a baryon. As a result, meson and baryon momen- 
tum distributions are proportional to the square and the 
cube, respectively, of the quark momentum distribution, 
This has led to the scaling of hadron elliptic flows with 
the valence quark number in the hadron [3, 0, Q , a novel 
scaling of heavy hadrons when the coalescing quarks have 
different masses and an enhanced baryon-to-meson 
ratio @, 1,0,0. 

However, the current formulation of the quark coales- 
cence model only applies to partons above a moderate 
transverse momentum scale pj_o, where the coalescence 
probability is small and most of these partons hadronize 
via other processes such as fragmentation. This limi- 
tation is due to the lack of the conservation of quark 
numbers of each flavor, or the lack of unitarity, during 
the coalescence process in the current quark coalescence 
formulation. For the bulk of partons in the quark-gluon 
plasma that are soft with transverse momenta well be- 
low the scale p±a, quark coalescence is expected to be 
the dominant hadronization process within the model, 
but a formulation including the quark number conserva- 
tion is necessary. Several earlier attempts to include uni- 
tarity have been made. In the algebraic coalescence re- 
hadronization (ALCOR) model, the quadratic and cubic 
scaling relations are used to describe meson and baryon 
multiplicities, respectively; and quark number conserva- 
tions are enforced by the normalization constants Q . An- 
other earlier study to include unitarity [Io| addressed the 



total multiplicities of different hadrons from quark coa- 
lescence and found that, as expected, hadron multiplicity 
is linearly proportional to the initial quark number. In 
a different approach, a multi-phase transport (AMPT) 
model [ll| converts nearby quarks into hadrons one-by- 
one and therefore satisfies the quark number conserva- 
tion; however, the model needs to be improved to better 
address the hadronization criterion that determines the 
average parton density at coalescence. 

In this study we incorporate the quark number conser- 
vation in the quark coalescence model so that the model 
can be applied to quarks at all transverse momenta. The 
general results as functions of quark transverse momen- 
tum are expected to yield the familiar quadratic and cu- 
bic scaling behavior [H, 0, [H, [|| for hard quarks. For 
soft quarks where the coalescence probability approaches 
one, the results should conserve the quark numbers of 
each flavor. 



2 The formulation including the quark number 
conservation 

Let us write f(p, x) = (2ir) 3 dN/(d 3 xd 3 p) as the phase 
space distribution of particles. A convenient formula for 
the coalescence production of meson M via a(3 — > M is 



E 



dN M (p) fda^p 
d 3 p 



:2n) 



A^A/I^p (q)| 2 /a(P«, X)f^p , x),(l) 



where p = p Q +P/3, q = p a — P/3, and gM is the statistical 
factor for forming the meson. The meson wavefunction is 
normalized as J c? 3 g| (q) | 2 = 1, and the first integration 
in Eq. (|TJ) runs over a 3-dimensional space-time hyper- 
surface. Assuming that the hadronization time scale is 
short so that the parton hypersurface during coalescence 
is unchanged and the hypersurface of formed hadrons 
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is the same as that of the coalescing quarks, we write 
E dN M (p)/d :i p=fd<j fJ 'p f j,fM(p,x)/(2ir) 3 , where the inte- 
gration is over the same hypersurface. The coalescence 
yield from a local hypersurface can then be written as 



/ M (p, x)=J d%g M |^ P (q)| fjpa,x)f^pp,x). (2) 

To include the quark number conservation, we need 
to take into account the fact that the quark number de- 
creases as quarks coalescence into hadrons. We thus con- 
sider the time evolution of the coalescence process. To 
differentiate from the time variable used to represent the 
hypersurface, let us use t c to represent the time variable 
of the coalescence process. For a local hypersurface at x, 
we write 

f'jjp, x, t c )=J d\C M {p, q, x, t c )f a (p a ,x, Qf^pp,x, t c ),(3) 

where f' M (p,x,t c ) = dfm(p, x, t c )/dt Cl and the phase 
space distributions such as f a {p a ,x,t c ) and fM{p,x,t c ) 
are time-dependent. The term Cm(p,<\, x,t c ) represents 
the probability for the coalescence production of meson 
M that could depend on x and time t c . 

When the internal momenta of the coalescing quarks 
in a hadron can be neglected, we can write 



C M (P, q, x, t c ) = c M (p, x, t c )S{q) 



(4) 



for meson M. When valence quarks have the same mass, 
Eq. reduces to 

/m(PM,I, tc) = C M {PM,X, <c)/a(Pq, X, t c )f^Pq, X, t c ),(5) 

where PM=Af 9 „P q with N qM =2, and c m (pm, x,t c ) repre- 
sents the coalescence coefficient for the meson. Similarly, 
for baryon productions via a + (3 + 7 — > B, we can write 

/s(PB,a;,tc)=CB(pB,a;,ic)/Q(Pq,^^c)//3(Pq,a;^ c )/7(Pq,a;,ic),(6) 

where PB=N qB Pq with N qB =3, and cb(pb,£,£ c ) repre- 
sents the coalescence coefficients for the baryon. 

The local conservations of quark numbers during co- 
alescence are given by a set of equations. An example 
is the following, which applies to quark flavor a when 
it appears as a constituent quark once in meson M and 
once in baryon B: 

/4(p qj x, t c )^NlJ' M (p M , x, t c ) - N 3 q J' B (p B ,x, t c ).(7) 

For a local coalescence process that starts at to and ends 
at tp, the initial conditions for coalescence are given by 

f a (Pq,X,t Q ) = /o(Pq,x), 

MPMii, to) = 0, /b(pb, x, t ) = (8) 

for quark flavor a, meson M and baryon B, respectively. 
In the following we mostly use the simplified notations, 
where the label x is omitted in all functions and the vari- 
able t r is rewritten as t. 



3 Consideration of one meson and one baryon 
species 

We consider the coalescence productions of one meson 
species and one baryon species via q + q — > M, q + q + 
q — > B, and q + q + q — > B. Let us assume zero baryon 
chemical potential for simplicity. For quark distributions 
at momentum p q , we have the following rate equations 

/m(PM,*) = CAl{pM,t)f q (Pq,t), 
f' B (.PB,t) =C B (pB,t)/ g 3 (Pq,0- (9) 

The conservation of quark numbers is given by 

4(p q , t)=-NlJ' M {p m ,t) - 3N 3 J' B (p B ,t). (fO) 
Therefore the quark distribution is given by 
/:(p q ,t)=-^ /C M(PM,i)/'(Pq,i)-3^s(PB,t)/„ 3 (p q ,i).(ll) 



3.1 Solutions when coalescence coefficients have the 
same time-dependence 

If the meson and baryon coalescence coefficients have 
the same time-dependence, we can define 



= W* B c B (p B ,t) 
{Pci) - Nl, CM(PM,t)' 



(12) 



Eq. (JTTJ) then has the following solution for / g (p q ,i) 

l + KPq)/o(Pq) \ ( fq(P^t) 



r(pq) In 
1 



/g(Pq, t) /o(p q ) 



1 +' , (Pq)/ 9 {Pq^)/ V /o(Pq) 

Im{pm, t), 



with Ij\/(pm, t) = Nq M / caKpm, 



it)du. 



(13) 



We first consider the limit r (p q ) — * 0, where meson 
productions dominate. We then obtain the solutions at 
the leading order in r(p q ) as 

^"I+MpSpm,*)' 



/a/{Pm ,t) : 



/o(Pg) J M(PM,t) 
Nl M [l + fdjP^)I M {pm,t)\ 



r (n ,y, /o(Pq) r (Pq)J 1 
MPB,t) ^ — jtto <1 



6N *b { [l+fdpMp^t)} 2 

Therefore the coalescence probability of quarks at mo- 
mentum p q is given by 



_ /_ x 1 /g(PqiM -1 l _ 

PwAPqj-l ^ l + f( ( Pq) I M ( PM ,t F ) 



(17) 



3 



If we define the scale p±o so that p coa ;.(Pq) = 1/2, p±o 
then corresponds to the quark transverse momentum 

where /o(Pq)-fM(PM,^) = 

For hard quarks, those with transverse momenta well 
above the scale p_lo, they mostly hadronize via other 
processes such as fragmentation, thus p CO ai. (Pq) <C 1, 
which is equivalent to /o(Pq)^Af(PM, t F ) *C 1 according 
to Eq. (fT7|) . This is consistent with the fact that /o(Pq) 
decreases rapidly as the transverse momentum increases. 
We then obtain the final hadron distributions from the 
coalescence of hard quarks as 

/mIpm,^) * f^)I M (PM,t F )/Nl r 



1M ' 



fB rd (PB,t F ) ~ / 3 (pq)/ B (p B ,tF)/(37V ( 3 
i — o at3 



with J B (p B , t) ee 3iV? / c s (p B , u)du 



(18) 



to 



Note that for hard quarks we have / 9 (p q ,£) — io(p q ), the 
time-integration of Eq. ([3]) then gives 



rtF 

/ c M (p,<i,x,t c )dt c ~ g M |* P (q)| s 
Jt 



(19) 



When the internal momenta of the coalescing quarks in 
a hadron are neglected, the above reduces to 



Similarly, we also have 



c B {p-B,t)dt ~ g_B 



(20) 



(21) 



for fast quarks. Therefore the functions -Zm(pm, t F ), 
-?b(Pb,£f) and r(pq) for fast quarks are constants that 
do not depend on cf>, the azimuthal angle of the momen- 
tum vector p q in the transverse plane of a heavy ion 
collision. We then have 



/M r W<F) « / 2 (Pq), 

f H B ar \p^t F ) OC / 3 ( P q). 



(22) 



These are just the scalingrelations of the previous quark 
coalescence model P, 0,M Hi- Consequently, when the 
quark momentum distribution has an azimuthal asym- 
metry that is dominated by the cos (20) term, the above 
scalings lead to the following quark-number scaling of the 
hadron elliptic flows (v 2 ) in ultra-relativistic heavy ions 
collisions hi: 



J 2,M 

har> 
J 2,B 



d (P±) 



2v 2 , q ( Pl jN qM ), 
3v 2>q (p±/N qB ). 



(23) 



We note that the fragmentation process, which domi- 
nates the particle yields at very high transverse momen- 
tum [f|, is not considered in this study. 



For soft quarks, those with transverse momenta well 
below the scale p±o, they mostly coalesce and thus 
PcoaiiPq) ^ 1, i.e., /o(p q )/M(PM,^F) > 1. The final 
hadron distributions from the coalescence of soft quarks 
are given by 



f B ° ft (PB,t F ) 



/ (Pq)/< f , 
/0(PqMPq)/(6^ 3 B ; 



(24) 



The above linear scaling of the dominant mesons is the 
result of including the quark number conservation in the 
formulation. Consequently, it leads to the following scal- 
ing of meson elliptic flows: 



soft, 

J 2,M 



\p±)^v 2<q (p±/N qM ). 



(25) 



On the other hand, the scaling of the non-dominant 
baryons is quadratic in this case. If r(pq) for soft quarks 
also does not depend on the angle </>, we have 



v 2°bXp±)- 2v %<i(p±I N <ib)- 



J 2,B 



(26) 



Comparing the scalings of Eqs. (|25H26f with those in 
Eq. (f!?3"|) , we find that the scalings of hadron elliptic flows 
with the quark elliptic flow for soft quarks are weaker 
than those for hard quarks. This can be understood 
from Eq. (|14|) . which says that the azimuthal asymme- 
try in the quark momentum distribution during the co- 
alescence process will be reduced from its initial magni- 
tude. For hard quarks this reduction is negligible because 
/o(p q )^/(PM, t) «C 1 is true throughout the coalescence 
process. For soft quarks, however, the azimuthal asym- 
metry in the quark momentum distribution will gradually 
vanish during coalescence, leading to weaker scalings of 
soft hadron elliptic flows with the quark elliptic flow. 

Secondly, let us consider the opposite limit r(pq)^oo 
by keeping cb(Pb,*) fixed while decreasing cj,/(pm,')- 
In this case, baryon productions dominate, and we obtain 
the solutions at the leading order in l/r(p q ) as 



/g(Pq,*): 
/a/(pm, t) - 
/b(pb,£): 



./o(Pq) 



v /l + 2/ ( ^(p q )/ B (pB,t)' 
ln[l + 2/ 2 (p q )/ j3 (p B ,0] 

2^V(Pq) 

,Mp q )( 1 i ) 



(27) 



The scale p±o in this case corresponds to the quark trans- 
verse momentum where /o(p q )/s(PB, ^f) = 3/2. For hard 
quarks well above the transverse momentum p±o, the fi- 
nal hadron distributions are the same as Eq. (fT8|) . How- 
ever, hadron distributions from the coalescence of soft 
quarks are given by 



f S M f \PM,t F ) 
f S B° ft (PB,t F ) 



ln[2/ 2 ( Pq )/ B ( PB ,<)] / [2AT 3 M r(p q )] 



/o(Pq)/(3iV 3 ). 



(28) 
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We see a linear scaling of the dominant baryons in this 
case due to the quark number conservation, but the scal- 
ing of the non-dominant mesons is different. Similar to 
the r(pq) — > case, from Eq. (|2"7|) we also see that the 
asymmetry in the quark momentum distribution of soft 
quarks will gradually vanish during the coalescence pro- 
cess, leading to weaker scalings of hadron elliptic flows 
with the quark elliptic flow for soft quarks than for hard 
quarks. 



3.2 Solutions when mesons dominate 

Without assuming that the baryon and meson coa- 
lescence coefficients have the same time-dependence, we 
now consider Eqs. ([9lllip in the limit that mesons domi- 
nate. In this case, the leading-order solutions of /q(p q ,i) 
and /m(pm,*) ar e still given by Eq. (fT4"|) and Eq. (fTS"]) . 
respectively. However, the baryon distribution is given 
by 



equations: 



f B {pB,t)c 



fo(Pg)cB(PB,u)du 

[l + /o(Pq)MPM,Ii)]' 



(29) 



For hard quarks well above the scale pj_o, the above is 
the same as the baryon distribution in Eq. (|18|) . For 
soft quarks, however, the scaling behavior of the non- 
dominant baryons depends on the coalescence dynam- 
ics, which are represented by the coalescence coefficients 
CAi(PM,t) and Cs(pB,i)- For example, if cb(pb,$ = 
e c A /(p M , *)/m(Pm,*) with e <C 1, we have 



/s(PB,t) 



e/o(Pq)^Af(PM,*) 



2Nl M [l+ffoJI M (p M ,t)Y 
For soft quarks, this leads to 

psoft. 



(30) 



/JTWij - e/^ Pq )/(2< f ), 
V2,Jp±/N qB ). 



soft/ \ 
V 2.B (P-U 



(31) 



On the other hand, if cb(pb,£) oc Cm(pm, t) as assumed 
in Eq. ()12p . the baryon distributions of Eq. fTS]) and 
Eq. (J2U) are reproduced, where soft baryons follow a dif- 
ferent scaling as given by Eq. ([26]) . Therefore the dy- 
namics of coalescence can be probed by studying hadrons 
from the coalescence of quarks below the scale p±Q. 



4 Charm mesons 

We consider the coalescence production of charm 
mesons via q + c — > D, q + c D, and c + c — ■> J/V'j m 
addition to the coalescence production of a light meson 
species via q + q — > M . Assuming zero baryon chem- 
ical potential for simplicity, we solve the following rate 



/m(Pm, t) = c M (p M , i)/ 9 2 (Pq> *)> 

/d(PD,*) = C£)(PD,i)/ 9 (Pq,t)/c(Pc,i), 
/V>(P-0I *) = C </>(P^ *)/c(Pc, *), 



(32) 



together with the following equations that represent the 
conservations of light quark and charm quark numbers: 



f q (p v t) = -<4(p M ,t) - f' D (pn,ty, 

f' c ( Pc ,t) = -f D (pn,t)/ : 



qD ■ 



'■c. D -N* M f^,t). (33) 



In the above, z qD — m q /(m c + m q ), z CD = m c /{m c + 
m q ), where m q and m c are the effective masses of the 
valence light quark and charm quark, respectively. The 
light quark momentum involved in the above equations 
is p q , and as a result the charm quark momentum in 
consideration is given by p c = Pqm c /m q 4]. We also 
have p D = Pq/z qr> and p^ = N qM p c . 

Assuming that light mesons dominate over D mesons 
in the light quark sector and D mesons dominate over 
J/ip mesons in the charm quark sector, we obtain the fol- 
lowing solutions at the leading order in / c o(Pc)//o(Pq) : 



Up, 



■J) ^ /co(p c )e Jt o , 



(34) 



where / c o(Pc) = /c(Po *o) represents the charm quark dis- 
tribution just before coalescence, and /q(Pq,i) is given by 
Eq. (fT4]) at the leading order. The light meson distri- 
bution /m(pm,£) is given by Eq. (fl~5|). and the D meson 
distribution is given by 



/ D (p D ,i)~/ c0 (p c )4fl-e J *o 



.(35) 



We now consider the case where the coalescence coef- 
ficients have the same time-dependence, and define 



^d(Pd) = 



7>(P</-) 



Cp(PD,t) 

c^lp^t)zl D Nl M 



C£)(pD,t) 

The solutions then simplify to 

/co(Pc) 



[l + /o(p q )/M(PM,i)] rD(PD) ' 



(36) 



(37) 



/D(PD,i)-/c^Pc)^<l- 



[l + /o(p q )/A/(pM,t)] r£,(PD) J ' 



U(P4,,t) 



/ e 2 o(Pc)r J D(p D )r^,(p^,) 
/o(p q ) [l-2r D (p D )]N*„ 

x {[l + ^Pq^^PM,*)] 1 ^^-!}. (38) 

According to Eq. (|37[) . the coalescence transverse mo- 
mentum scale for charm quarks is given by p'j_ m c /m q , 
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where p' ±0 corresponds to the light quark transverse mo- 
mentum where /o(Pq)^/(PM, t) = 2 1 / r °^ - 1. Thus 
p' ±Q could be quite different from p±o, the coalescence 
transverse momentum scale for light quarks. 
For the coalescence of hard quarks, we have 

/D ar< W^ io(Pq) /co(Pc) zl D r D (pr>) I M (pm, t F ) , 
4 Qr tp^tF)^/ c 2 (Pc)r-v,(p^)r D (p D )/M(p M3 tF)/^ 3 M j(39) 

which reproduce the previous scaling relations jij. For 
soft quarks, we first note that, according to Eq. ([57)) and 
the solution of fqiPqJ) in Eq. (TH)) . the assumptions that 
we made in obtaining the solutions of Eqs. (|37H38p are 
valid when ?\d(pd) > 1- If ?\d(pd) < 1, the charm quark 
distribution decreases with time at a slower rate than 
the light quark distribution, and therefore the number 
of charm quarks would be comparable with that of light 
quarks after a certain time. However, this cannot happen 
as long as [/c/pJ/mCpm,^)] 1 ^^ < /o(p q )//co(Pc), 
and we limit our discussion to this case. For soft quarks 
with /a(Pq)-fM(PM,*F)»l, we then have 

fD ft (Pr>,t F ) - /co(Pc)*c D , 
.soft, v / c 2 (Pc)^(p^)r g (p D ) 

^^-/W^tPo)-!]^/ " - 1 ' 

/co(Pc)^(p^)^(Pd) , * . . , /0 

* /o( Pq )[2,.( PD )-l]^ fOT 1 ^ >1/2 ' 



/ c 2 o(Pc)^(p )r D (p D ) 

77- w r3 m[/n(Pq)iM(PM,lF)J 

JdPl)^q M 

for r D =l/2, 
/co(Pc)^(P^) r s(PD) 



/^ Pq )[l-2r D (p D )]iV3 M 
x [Zo(p q )/M(pM,tF)] 1 " 2rD(PD) for r I? <l/2.(40) 

Note that the above solution of fl° (p^tp) for 1 > 
ro>l/2 is valid if t\d is not very close to 1/2 so that 
LMp q )//co(Pc)] (2rD ~ 1)/(1 - r£,) » 1 is satisfied, and the 
above solution for ro<l/2 is valid if rp is not very close 
to 1/2 so that [/o(p q )^Af(PM,iF)] 1 ~ 2rD » 1 is satisfied. 

We find that, contrary to naive expectations, the dis- 
tribution of soft D mesons only depends on the charm 
quark distribution. As a result, we have 



V2 dXp±)-V2A Z c D P±), 



(41) 



which is different from the result for hard D mesons 



(42) 



Furthermore, the distributions of soft J/ip mesons 
are proportional to / 2 o(Pc)//o(Pq) with a being positive 
(modulo a logarithmic dependence on /o(p q )), and this 
leads to 

v s 2 %%P±)^2v 2 AP±/ N qM) - av 2 ,q{p±m q /Nq M /m^.(A3) 



For example, if there is no charm elliptic flow, the elliptic 
flow of J/ip mesons would be negative since light quarks 
have a positive elliptic flow. We may understand this cor- 
relation from Eq. (137p , which shows that the charm quark 
distribution during the coalescence process is inversely 
correlated to the light quark distribution /o(p q )- If there 
is no azimuthal asymmetry in the charm quark momen- 
tum distribution from the parton phase, fdj>c,to) does 
not depend on the azimuthal angle of the charm quark 
momentum p c . After time to-, however, Eq. (|37|) shows 
that the charm quark momentum distribution f c (p c ,t) 
will have an azimuthal asymmetry that has an opposite 
sign to the asymmetry in the light quark momentum dis- 
tribution. Because the formation rate of the J/ip meson 
distribution at any given time during coalescence is pro- 
portional to the square of the charm meson distribution 
at that time, J/ip mesons thus have an elliptic flow that 
has an opposite sign as the light quark elliptic flow. 



5 Discussion 

The final hadron distributions that we derive in this 
study are the distributions just after quark coalescence. 
Later rescatterings in the hadron phase, which could sig- 
nificantly affect hadron momentum spectra and particle 
ratios, have been neglected. The distributions are only 
derived for a local hypersurface, therefore the effect from 
integrating over the full hypersurface is not studied. Nu- 
merical estimates of quantities such as the initial quark 
distributions and coalescence coefficients need to be done 
to obtain the relevant solutions as well as the coales- 
cence momentum scales for relativistic heavy ion colli- 
sions. Also missing is the effect of the internal momenta 
of coales cing quarks, which may be important for soft 
hadrons [131 ] - Furthermore, resonance decays have been 
shown to affect the elliptic flow, especially that of pions 
[l3j : however, we have not considered multiple hadron 
species including resonances. 



6 Conclusions 

We have developed a quark coalescence model that 
includes the quark number conservation. We find that 
scalings between hadron and quark momentum distribu- 
tions depend on the quark transverse momentum. For 
hard quarks, our general results reproduce the usual scal- 
ing relations. For softer quarks, however, hadrons that 
dominate the quark number of each flavor exhibit a lin- 
ear scaling due to the quark number conservation, while 
the scalings of non-dominant hadrons are sensitive to the 
coalescence dynamics. Charm mesons are also studied. 
Contrary to naive expectations, we find that the distribu- 
tion of soft D mesons does not depend on the light quark 
distribution, while the distribution of soft J/ip mesons is 
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inversely correlated to the light quark distribution. These 
quark-hadron scalings as functions of particle momentum 
can be used to test the quark coalescence model as an ef- 
fective hadronization model for partonic matter, and the 
confirmation would provide new evidence for the forma- 
tion of the quark-gluon plasma in relativistic heavy ion 
collisions. 
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